Complete analytic solution of the geodesic equation 
in Schwarzschild—(anti) de Sitter space—times 
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The complete set of analytic solutions of the geodesic equation in a Schwarzschild-(anti) de Sitter 
space-time is presented. The solutions are derived from the Jacobi inversion problem restricted 
to the theta-divisor. In its final form the solutions can be expressed in terms of derivatives of 
Kleinian sigma functions. The solutions are completely classified by the structure of the zeros of the 
characteristic polynomial which depends on the energy, angular momentum, and the cosmological 
constant. 
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Introduction and motivation The Einstein equations 
with a cosmological constant 

today are taken as very successful model to incorporate 
the observation of an accelerated expanding universe [l| . 
Within this theory the most simple solution of an isolated 
mass is the spherically symmetric Schwarzschild-de Sit¬ 
ter solution (e.g. i) 

ds^ = adt^ — a~^dr^ — r^{d0^ + sin^ Odip ), (1) 


where a = 1 — ^ (we take c = G = 1). This 

Schwarzschild-de Sitter metric is characterized by the 
Schwarzschild-radius rs = 2M related to the mass M 
of the gravitating body, and the cosmological constant A 
which today is assumed to have the value ^ 10“®^ m“^ 
1|. For a general discussion of this metric, see e.g. 
2, Q . This may serve as a model for the Solar system 
in an expanding universe. Recently such a model has 
attracted some interest in order to discuss the possibil¬ 
ity to observe the influence of A on moving masses on 
the scale of the Solar system Q. This requires a discus¬ 
sion of the geodesic equation and its solutions in such a 
Schwarzschild-de Sitter space-time. Exact solutions of 
the geodesic equation are also of interest in the calcula¬ 
tion of gravitational wave templates for Extreme Mass 
Ratio Inspirals (EMRI) [^. There is further interest to 
understand explicitly the structure of geodesics in the 
background of black holes in anti-de Sitter space in the 
context of string theory and the AdS/CFT correspon¬ 
dence. In addition, recently there also has been a lot 
of work dealing with geodesics and integrability in black 
hole backgrounds in higher dimensions, also in the pres¬ 
ence of a cosmological constant @ . 

Here we present explicitly the full analytic solution of 
the geodesic equation in a Schwarzschild-(anti)de Sit¬ 
ter space-time and their complete classification. So far, 
explicit solutions for the geodesic equations are only 
known for the Schwarzschild solution [Tj, the Reissner- 
Nordstrom solution H , and for the Kerr and Kerr- 
Newman space-time [8|. For a specialized case orbits 
in a Schwarzschild-(anti) de Sitter space-time have been 
presented in [^. 


The geodesic equation We consider the geodesic equa¬ 
tion 


r /i 1 dxP dx'^ 

dg2 ’ 


( 2 ) 


where ds^ = g^^dx^dx’' is the proper time and {pa} is 
the Christoffel symbol in a space-time given by the met¬ 
ric O- The geodesic equation has to be supplemented 
by the normalization condition = e where for 

massive particles e = 1 and for light e = 0. 

Due to the spherical symmetry we can restrict our 
consideration to the equatorial plane. Furthermore, due 
to the conserved energy E and angular momentum L 
the geodesic equation reduces to one ordinary differen¬ 
tial equation 


1 /dr'. r / N 

b : 


2 \ds 

with the effective potential 

K.ff(r) - - (^--AL -e- + —- 

A substitution u = r^/r gives 
2 


( 3 ) 


-AE 


( 4 ) 


du 

dip 


1 


1 


= — u^ eXu + — e) + -p+ ^Xp— (5) 

with the dimensionless parameters 


A := -^ > 0 , p := >0 and p := Ari . (6) 

We rewrite (O as 

with 


P5{u) := 


u‘^+eXu^+ 




( 8 ) 
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For light e = 0 the situation simplihes considerably. 
In this case the solutions can be given in terms of the 
Weierstrass p function 0 


r{p) 


_ 


(9) 


as well as the period matrices ( 2 a;, 2 a;') and (27^, 2rj') 




(15) 


with the Weierstrass invariants 


Here {ai,bi \i = 1,2} G Hi{X,'L) is the homology basis 
of closed paths where integration does not evaluate to 
zero. Finally we introduce the normalized holomorphic 
differentials 


where uq = rg/b and b is the distance of closest approach. 

Now we discuss e = 1. Separation of variables in d?]) 
yields the hyperelliptic integral 


ip-ipo 



u'du' 

/W) ’ 


( 11 ) 


where uq = u{(po). We explicitly solve this integral by 
using a method based on the inversion problem of hyper- 
ellip tic Abelian integrals [Il| . We use this general ansatz 
[l2l | in order to obtain an explicit one-parameter solution 
of the two-parameter inversion problem by restricting the 
problem to the so-called theta-divisor. This procedure 
has been applied previously to the problem of the double 
pendulum [l3|. The resulting orbits can be completely 
classified. 

In solving (jllD there are two major issues which have 
to be addressed. First, owing to the two branches of 
the square root the integrand is not well defined in the 
complex plane. Second, the solution u{ip) should not 
depend on the integration path. That means that if 7 
denotes some closed integration path and 


udu 


= bJ 


( 12 ) 


then 


dv := ( 2 a;) ^dz , 


dz = 


dzi 
dZ2 


(16) 


The period matrix of these differentials is given by 
( 12 ,t), where I 2 is the 2 x 2 unit matrix and r is de¬ 
fined by 


T := oj . 


(17) 


It can be shown 1^ that this normalized matrix r always 
is a Riemannian matrix, that is, r is symmetric and 
is positive definite. 

The Riemannian surface of = P 5 {x) has 4 inde¬ 
pendent closed paths, each corresponding to a period of 
the functions dehned on these surfaces and, hence, to a 
period of the solution u of m- In order to construct 
functions with four periods, we need the theta function 
7 ? : C2 ^ C, 


7?(z; t) ■= 


^ ^ (Tm-|- 2 z) 

rh^l? 


(18) 


which is a holomorphic function on (D^. The theta func¬ 
tion is already periodic with respect to the columns of I 2 
and quasi-periodic with respect to the columns of t, i.e. 
for any n € it holds 


-d(z + l 2 n;T) =-d(z;T), (19) 

i9(z + rn; r) = r). ( 20 ) 


(p-<po -UJ = 



u'du' 


(13) 


should be valid, too. Therefore u has to be a periodic 
functions with period oj. 

The analytic solution These problems can be solved 
if (HU) is considered on the Riemannian surface of = 
P 5 {x). For the construction of periodic functions on this 
Riemannian surface we introduce the canonical basis of 
the space of holomorphic differentials dzi and of associ¬ 
ated meromorphic differentials dvi by 


dzi 

dri 


dx 


VMx) 

3x^ — 2x^ -I- Aa; 


dz2 ■■= 


dx, dro := 


xdx 

y/P^ix) ’ 
x^dx 

WP^ix) ’ 


(14) 


The solution of m can be explicitly formulated in 
terms of sigma functions ct : C, 

a{z) = (^{2 uj)-^z + + > 

( 21 ) 

where the constant C can be given explicitly 1^ but 
does not matter here. Based on this a are the generalized 
Weierstrass functions 

dd a,{z)aj{z) - a{z)a,j{z) 

ay ^’ 

( 22 ) 

where Ui denotes the derivative of the sigma function 
with respect to the i—th component. 

We need the 2-dimensional Abel map 


: S'^X Jac(A), 


i-t 


dz+ / dz (23) 
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FIG. 1. The zeros of P^iu) in a (/r, A)-diagram for different values for the cosmological constant is along the a;-axis, A along 
the j/-axis). The grayscales encode the numbers of positive real zeros of the polynomial P 5 : black = 4, dark gray = 3, gray = 
2, light gray = 1, white = 0. For A = 0 these zeros have been extensively discussed in [J. 


from the 2nd symmetric power of the Riemannian surface 
X to the Jacobian Jac(Jf) = C^/F of X, where F = {wv-\- 
uj'v'\v,v' e I?} is the lattice of periods of the differential 
dz. Jacobi’s inversion problem is now to determine x for 
given ip from the equation ip = Ax„(x), i.e. 


(^1 = 

= 


dz 


xo 

xi 


zdz 


L 


dz 


Xq 

Xl 


zdz 






(24) 


We want to restrict now Jacobi’s inversion problem to 
the set of zeros of i?(- + where is the vec¬ 

tor of Riemann constants, see This set is called a 
theta-divisor and denoted by ©5; . This restriction is 

necessary since in our case A~^ = ip should be the (one¬ 
dimensional) physical angle and the theta-divisor is a 
one-dimensional subset of S'^X. 

The solution of Jacobi’s inversion problem (IMl) can be 
given in terms of generalized Weierstrass functions. The 
components of the solution vector x = {xi,X 2 y are 


Xi+X2= 4:p22(ip) , 
Xl X2 = -4pi2(<p) ■ 


(25) 


Now we reformulate (l24ll to c/) = ^00 (x) with (j) = 0 — 
dz and restrict the problem to the theta-divisor by 
the limiting process 


Xl = lim 


X 1 X 2 


X2^oo Xl + X2 

fj(‘Poo)(^12(0 00 ) - CXD )O'2(0 00 ) 

ct|(<Poo) - Cr((^oo)a’22(<Poo) 


(26) 


where 


(^00 = lim 0 = ^oo(^oo), ^00 = (a^i, 00) 

ai2—^00 


(27) 


It can be shown, see [1^, that (2a;) ^Aoo(xoo) is an el¬ 
ement of the theta-divisor 0^ , where Koc = I'j + 


[HI- Thus, (7{0ao) = 0 and there is a functional 
relation between poo,! and fao, 2 - The solution to the 
geodesic equation in Schwarzschild-de Sitter space-time 
is then given by 


, ^ rs rs ( 72 (^ 0 ) 

rm = —rw = — = -rs— 

u{ip) Xl cri(ipQ) 


where 00 = (^ ) . Here ip'^ is given by 

poo 

^'0 = ^ 0 + / dz2 

d Un 


(28) 


(29) 


and Pi is chosen in such a way that [2uj)~^0q is an ele¬ 
ment of the theta-divisor ■ 

The classification of the solutions Having an analyti¬ 
cal solution at hand we can explore the set of all possible 
solutions in a systematic manner. The shape of an orbit 
depends on L, and A or, equivalently, on A, /x, and p 
(rs can been absorbed through a rescaling of the radial 
coordinate). Since r should be real and positive it is clear 
from (|4|) that physically acceptable regions are given by 
those u for which E > Ves. The zeros of P5 are related 
to the points of intersection of E and I4ff, and a real and 
positive P5 is equivalent to if > 14®. Hence, the number 
of positive real zeros of P5 characterizes the form of the 
resulting orbit. A detailed discussion is presented in 0. 

If we denote by ei,...,e„ the positive real zeros, 
then it follows that the physically acceptable regions are 
given by [0, ei], [e2, 63], ..., [e„, 00] if n is even and by 
[ei, 62 ],..., [e„, 00] if n is odd. With respect to r we have 
the following classes of orbits: 


(i) the region [0, ei] corresponds to escape orbits, 

(ii) the region [e„,oo] corresponds to orbits falling into 
the singularity, i.e. to terminating orbits, and 

(iii) the regions [ci^Ci+i] correspond to bound orbits. 


The case of no positive real zero corresponds to a particle 
coming from infinity and falling to the singularity. 

The number of real zeros are shown in Fig. [T] A non¬ 
vanishing cosmological constant dramatically changes the 
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(b) /i = 1.1, A = 0.2, 
A = -10-5km-2^ 
ro = 185.37km 


FIG. 2. Orbits with no analogue for A = 0. (a): Reflection 
at the A-barrier (b): bound orbit with perihelion shift which 
would be an escape orbit for A = 0. The scale is in km. 


possible structure of orbits. While for negative cosmo¬ 
logical constant there are only bound orbits or orbits 
terminating in the singularity, for vanishing A we have 
terminating, bound as well as escape orbits. For A > 0 
we have in addition orbits which will be “reflected” at 
the potential barrier related to positive A. In Fig. [^two 
orbits without analogue for A = 0 are shown. 

For bound orbits it is possible to directly identify the 
perihelion shift Aperihei by the difference of the period¬ 
icity of the solution r{ip) from 27r. Let us assume that 
the bound orbit corresponds to the interval [cfc, ek+i] and 
that the path ai surrounds this interval. Then 

r^^+1 xdx 

Aperihei = 27r — 2u;2i = 2 / , . . ■ (30) 

•fe/c 


Summary and outlook We derived the complete set of 
solutions for the geodesic equation in a Schwarzschild- 
(anti)de Sitter space-time. A more detailed discussion 
of the approach will be published elsewhere 0. Beside 
the importance of such solutions mentioned in the in¬ 
troduction this result may also be applied to the study 
of the Pioneer anomaly [l^ which has been speculated 
to be due to the influence of the cosmological expan¬ 
sion. Although the cosmological constant is very small, 
its influence on nearly critical orbits is much stronger. 
As the anomaly is also small, a definite quantitative 
answer to the question whether the cosmological con¬ 
stant has an influence on the Pioneer orbits requires 
an analytical solution. Our result may also be applied 
to the motion of stars in the outer regions of galax¬ 
ies. Furthermore, the above techniques apply to the 
geodesic equations in higher dimensional Schwarzschild, 
Schwarzschild-(anti)de Sitter, Reissner-Nordstrom and 
Reissner-Nordstr6m-(anti)de Sitter space-times [ 13 . The 
use of this method in more general space-times, namely 
the Plebanski-Demianski class without acceleration is 
under way 17|. In these cases polynomials of 6th order 


appear which can be treated with the above method. 

We note that our method is not capable to solve equa¬ 
tions of motions based on a polynomial of 7th or higher 
order which appear, e.g., in the geodesic equation in 
Schwarzschild space-times in higher than 7 dimensions. 
This will entail an Abel map with a preimage of dimen¬ 
sion three or higher, which could only be reduced by one 
dimension with the presented method. 
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J. Kunz, O. Lechtenfeld, D. Lorek and in particular P.H. 
Richter for fruitful discussions. Financial support of the 
German Aerospace Center DLR and of the German Re¬ 
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